Direct Counting Analysis on Network Generated by Discrete Dynamics 
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A detail study on the In-degree Distribution (ID) of Cellular Automata is obtained by exact 
enumeration. The results indicate large deviation from multiscaling and classification according 
to ID are discussed. We further augment the transfer matrix as such the distributions for more 
complicated rules are obtained. Dependence of In-degree Distribution on the lattice size have also 
been found for some rules including R50 and R77. 
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Discrete dynamical system in general can generate very 
complicated dynamics even though the rule or the equa- 
tion of evolution for the system can be apparently simple. 
Examples include the synchronization of Pulse-coupled 
Oscillators, chaos arises from discrete mapping and more 
interestingly the simple Cellular Automata (CA) intro- 
duced by von Neumann almost half a century ago [l|. 
One of the important problems in discrete dynamical sys- 
tem is the classification of it's general behavior. Gener- 
ally speaking this is a hard question to be answered since 
the classification requests a full understanding of the dy- 
namics itself. Even for simple system such as CA which 
evolves according to simple rules, a complete classifica- 
tion is still controversial since the work of StevenWolfram 
in the 80' of last century 0. Recently, A. Shreim et al. 
have approached the CA problem by employing concepts 
originated from network analysis Q. They classified CA 
by analyzing the local and global properties of the net- 
work generated by the dynamics of CA. Their approach 
certainly provides new way to address this interesting 
and difficult problem. In this work we follow up their 
approach and investigate further on this problem. Here 
we would only concentrate on the local analysis and will 
return on the global property later. 

For any discrete dynamical system, the configuration 
space can be represented by discrete nodes. The time 
evolution of the system in each time step is presented by 
a connected link starting from one node to its dynamical 
successor. As a result, all the trajectories of the dynam- 
ical system generated a directed network. Since CA are 
deterministic systems, each node possess one single out- 
going link. However, due to the existence of fixed points 
and periodic solutions, the dynamics of CA in general is 
irreversible and the number of preimages of a state can 
be larger than 1. The in-degree (ID) of a node is de- 
fined as it's number of preimages and obviously in CA 
the in-degree is a local property. In [3[ the authors have 
performed analysis on the ID distribution for various one 
dimensional CA both analytically and numerically. Due 
to the simplicity of Rule 4 (R4), by using a transfer ma- 
trix (T) approach together with combinatory analysis on 
preimages, they obtained an approximation solution of 
the probability distribution P(k) with k denoting the in- 
degree. Their results were presented on the log — log plot 
as curves instead of straight lines and seem to suggest 



the existence of multiscaling distribution. In our work, 
instead of using approximation, we re-analyze this prob- 
lem by performing direct counting of the in-degree for 
various rules with different lattice size L. What we have 
found is that the distribution is not as simple as what 
have been obtained in [3J. The plan of this report is as 
follows. A brief summarization on nomenclatures of CA 
and network are given and the results on ID of [|[ are 
briefly discussed in section [U The direct counting results 
of various rules are reported in sectionllTl where some fur- 
ther studies on the ID are also treated. The extension of 
transfer matrix for more accurate analysis on other rules 
and its implication are treated in Section Mil A brief 
summary of our findings is provided in the final section. 



I. INTRODUCTION 

A one-dimensional Cellular Automaton is defined on 
a spatial lattice of L sites where L can either be finite 
or infinite. In this work we will keep L finite and the 
state of each site is in one of the g states at any time 
t. Each site follows the same prescribed rules for up- 
dating. For elementary CA, the number of neighbor- 
hood of each site is 2. The state of CA starts out with 
arbitrary initial configuration which is represented by 
'S'(O) = { s 0i si, S2, • • • , sl} where Si can be in any one 
of the g states. The configuration of system at time t 
is denoted by S(t) — {sq, s\ , s|, . . . , s l L }. In this work 
the state of each site is restricted to g = 2 such that 
Si € {0, 1}. Let R be any one of the CA, the value of 
Si at t + 1 time step is set equal to #(s-_ 1 s-s- +1 ) which 
is equal to or 1 according to the rule of interest. Ap- 
parently for elementary CA there are 256 rules. Follow- 
ing Wolfram's notation each rule can be assigned with a 
number given by: 



#(000) + 2#(001) 



-2 2 #(010) 
-2 5 #(101) - 



-2 3 #(011) 
2 6 #(110) 



- 2 4 #(100) 
2 7 #(111). 



For example, Rule 237 (R237) is the one with #(100) = 
0, #(001) = and all others R(sis 2 s 3 ) = 1. Based 
on large number of numerical studies, Wolfram has sug- 
gested that CA can be classified into four classes [2] . Cel- 
lular Automata within each class have the same qualita- 
tive behavior. Starting from almost all initial conditions, 
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trajectories of CA become concentrated onto attractors, 
the four classes can then be characterized by their at- 
tractors. According to Wolfram's classification, classes I, 
II and III are roughly corresponding to the limit points, 
limit cycles and chaotic attractors in continuous dynam- 
ical systems respectively. More precisely their respective 
long time limits are: (I) spatially homogeneous state, 
(II) fixed (steady) or periodic structure and (III) chaotic 
pattern throughout space. The fourth class of CA be- 
haves in a much more complicated manner and was con- 
jectured by Wolfram as capable for performing universal 
computation. For finite one-dimensional CA, the lattice 
is arranged on a circle with periodic boundary conditions. 
Such Cellular Automata have a finite number of states 
N = 2 L , and as a result, after sufficient long time evo- 
lution the system must enter the state which can either 
be homogenous, steady or periodic state. Therefore, the 
class III and IV attractors do not exist in finite CA. How- 
ever, it is also known that classes III and IV are hard to 
distinguished in some cases and resulting with disputed 
classification. In some sense this is due to the lack of a 
sensible way for defining complexity which is the crite- 
rion for class IV. In Q A. Shreim et al. studied networks 
generated by CA, they claimed to have found highly het- 
erogeneous state space networks for classes III and IV in 
contrast to the networks generated by classes I and II. 
One of the characteristic of state space is the in-degree 
distribution P(k) which is defined as Q 



P(k) 



S 
N' 



(1) 



where L is the size of the lattice, N = 2 L is the total num- 
ber of state and £ is the number of state with in-degree k. 
However in their analysis on the local properties, approx- 
imation has been applied in the analytical calculation of 
ID which shows nontrivial scaling. For completeness we 
briefly summarize their method of calculation in this sec- 
tion. They introduced the transfer matrix T which maps 
each pair Sj_iSj onto the pair SiSi+i, and such T^ s ' could 
be used to characterize all the preimages of hubstates (A 
hubstate is the state which has the maximum in-degree 
k m .) The rows and columns are order as "00", "01", "10" 
and "11" accordingly. For example the T matrix of R4 
is given by: 



110 
1 
110 
11 



(2) 



It is easy to check that the maximum in-degree k m is 
related to T as 



kn 



Tr{T L ). 



(3) 



The scaling of k m with the lattice size can be obtained by 
using the largest eigenvalue A m of T, namely, k m m N v 
with v — log 2 A m . Furthermore, for R4 we assigned n 



as the number of isolated l's, for example two config- 
urations (010000) and (010100) have n = 1 and 2 re- 
spectively Then by assuming a one-one corresponding 
relation between k and n, the following expression was 
proposed for the in-degree distribution: 



P(k)dk 



Q,(n)dn 
2 L ' 



(4) 



A multiscaling result is then obtained and given by the 
following expression Q: 



V = -1 -x + log 2 



(I-,) 1 



e £ (l -2c) 1 - 



2c 



(5) 



where y = logP(fc)/log N, x = logfc/logiV and e = n/L. 
The result of Eq. §5$ was ploted in @ and are apparent 
curvature appeared on the curve as such a conclusion of 
nultiscaling method than finite-size effect was proposed. 
Similar results of multiscaling in other rules were also 
reported in their article by using the same approxima- 
tion approach. To compare the results of Q with exact 
calculation we have performed the exact enumeration of 
R4 with different L. For the lack of space, in Table U 
we only list all the in-degree's for L — 12. The same 
characteristic also exist for other values of L. One no- 
tices that for the same k there corresponds more than 
one value of n. For example, for k = 12, n can either be 
2 or 3. Therefore the meaning of dn/dk is ambiguous. 
One might argue that the multi-valueness of n is an arti- 
fact of small L. However we have checked with L = 10 4 
where different states with the same k were found. In 
fact, for k = 1, states of different n can be constructed 
easily. For example, for k = 1, we have found at least 
2 states (000100010001 . . .) and (010101 . . .) which cor- 
respond to n = 2500 and 5000 respectively. Since the 
results were obtained for small k, there is a possibility 
that dn/dk might be well-defined for large k. To refute 
this reasoning, we have performed the large k analysis 
for different L. For L = 12, k m = 852 with n = and 
the second maximum in-degree ki m = 114 with n = 1; 
even in L = 25, we got k m = 1276941 with n = and 
&2m = 170625 with n = 1. The results indicate that for 
Art = 1 the Afc = k m — k2m is not a small quantity. 
In fact Ak rises exponentially as L getting larger. As a 
result, dk can never be infinitesimal implying the mean- 
ingless of dn/dk. Thus the simple multiscaling results in 
[3| is in doubt. 

The result of Eq. @ is presented by the solid line in 
Figure [1] By substituting the corresponding k m in Eq. 
([5]), one has y — —1.8. For R4 with L = 25, the results 
of n — k relation are provided in Figure [TJ There is only 
one hubstate with n = in such system, consequently 
the exact P(k m ) equals to 1/2 L Q which implies y = — 1 
independent of L instead of —1.8 as obtained by the ap- 
proximation calculation. Obviously the approximation 
for getting Eq. ([5]) is inadequate. In fact the above dis- 
crepancies prompted us to re-analyze this problem. 
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TABLE I. The relation of in-degree k and n for R4 with 
L = 12. 



TABLE II. The classification for different rules with A„ 
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FIG. 1. The relation of n and in-degree k for Rule 4. 
solid curve represents Eq. (6) in 
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II. NETWORK CHARACTERISTICS BY 
DIRECT ENUMERATION 



As suggested in @ , local properties such as the scaling 
of maximum in-degree with lattice size and the in-degree 
distribution might be useful for characterizing networks 
and hopefully also provide some help to classify CA. How- 
ever, from previous discussions it is clear that approxi- 
mation calculations can lead to inadequate conclusions 
such as multiscaling of ID being an example. In order to 
have a better understanding on this kind of network, it is 
unavoidable to address the problem with exact enumera- 
tion which will be the approach adopted in this work. We 
have studied more rules from different classes suggested 
by Wolfram and the in-degree distributions are eval- 
uated exaclty by direct counting method. The transfer 
matrix approach proposed in [3| has been used in this 
work extensively. 

It is important to note at the outset that not all rules 
can have a transfer matrix representation. The following 
matrix represents possible candidate of T for elementary 



An 



Rules 



1.32472 
1.46557 
1.61803 
1.75488 
2 



0.4056 

0.551 

0.6942 

0.811 

1 



110,62 

232,33,6 

36,18,5 

4,32 

0,255 



CA: 



a 2 

«3 B4 

a 5 a 6 

07 a§ 



(6) 



where a^, i = 1 ~ 8, can take on either 1 or 0. In fact 
there are totally 256 different matrices represented by T, 
however it is also known that not all of them can rep- 
resent the 256 elementary CA. This is partly due to the 
fact that not all rules can have a T matrix representa- 
tion. In particular, by diagonalizing all possible T, we 
have obtained only 9 different A m : 1, 1.22074, 1.32472, 
1.38028, 1.46557, 1.61803, 1.75488, 1.83929 and 2. Some 
of the rules corresponding to different A m are listed in 
Table El 

In what follows, all ID are evaluated by direct count- 
ing. To set the stage we present the in-degree distri- 
butions of various rules with different L in Fig. [5J the 
results of R4, R22 and R110 are plotted. The results of 
each rule seem to show the same distribution for three 
different values of L. Especially for small k the distri- 
butions are in good agreement for different L with all 
the points fall on a simple curve with non- vanishing cur- 
vature. As k increases, the distribution becomes more 
spread out but still maintains the same distribution for 
various L. These observation seems to suggest that the 
ID of this kind of network is scale independent. More- 
over, since these rules belongs to different classes and the 
scale independent of the distributions seem to be univer- 
sal for all four classes of CA. Similar results can also be 
obtained for other rules and will not be shown explicitly. 
Even though the L independence of ID is very sugges- 
tive, however a definite answer would require a full scale 
analysis with larger lattice size. Work along this line is 
now being studied and will be reported separately. On 
a different front, with knowing all the possible values of 
A m , it is natural to wonder if rules with the same A m 
are correlated. To show whether such correlation exists 
or not, results of the same v = log 2 A m are plotted in 
Fig. [3] and EO For v = 0.4056, the ID of R193, R147, 
R110, R137, R124 and R54 are given in Fig. [3J all of 
which show similar distribution. Interestingly these rules 
belong to the same class, namely class four according to 
Wolfram's scheme. For clarity, in Fig. [4] we have also 
plotted the ID of other rules of class IV. They also share 
the same feature as given in Fig. [3] Similarly, the re- 
sult for other are given in Fig. Ufa) with v = 0.551, the 
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corresponding rules are all of class II [5J and the ID of 
each rule also shows a universal distribution. However 
such good correlation between ID and v is not conclu- 
sive. Furthermore rules from different classes can have 
the same v. For example, the plot given in Fig. [5jb) for 
the in-degree distributions of various rules corresponding 
to v = 0.6942 shows diverse patterns. It is noted that 
in this plot, R36 and R5 are class II, whereas R18 and 
R126 are classified as class III Q. To make things worse, 
the plot for various rules of class II is given in Fig. EJc), 
where apparently the patterns get more diverse. More 
importantly, those rules belong to different z/'s. As a re- 
sult, one might conclude that there are good correlation 
for class IV which possess v — 0.4056, apart from that 
the ID of other classes do not share a universal pattern 
and also the value of v can not be used as a indicator for 
classification. Thus the local properties such as v and 
in-degree distribution are not useful for CA classification 
as suggested in [3] . 



III. AUGMENTATION OF T MATRIX 

The T matrix approach for calculating maximum in- 
degree is quite efficient and there are also cases in which 
the in-degree can also be evaluated exactly. However, as 
clearly stated in the previous section that 4 x 4 T matrix 
do not always exist for elementary CA and therefore it is 
interesting to see how modification on the transfer natrix 
scheme can be extended to other rules. In this section an 
extension is suggested by augmenting the 4x4 matrix to 
a 8 x 8 matrix. Instead of mapping each pair s\_ x s\ onto 
the pair s-_ 1 s-, T will map s*_ 2 s*_ 1 s| into the triplet 

s t i _ 1 s t i s t i+1 and hence define a 8 x 8 T matrix which could 
be used to describe the preimage more efficiently. The 
basis states are 

1 1) e {(in), (no), (ioi), (ioo), (on), (oio), (ooi), (ooo)}. 

Fig. [5] shows the definition of 8 x 8 T matrix. 

There are some merits in this approach. First of all, 
for elementary CA where all the rules of evolution involve 
three adjacent cells, obviously such extension certainly 
respect the evolution dynamics and hence include more 
information of the dynamical system. Secondly, the aug- 
mentation of T indeed resolves some of the cases which 
can not be done with 4x4 matrices. Moreover, this 
approach can also be applied to other more complicate 
CA including the model with next nearest-neighbor in- 
teraction, we will illustrate that by an example in this 
section. 

It is easy to see that whenever the 4x4 transfer matrix 
can be defined then there is always a corresponding 8x8 
matrix however the reverse is not always true. For the 
case when both descriptions exist then a natural ques- 
tion about the uniqueness of the results arises, namely, 
can both 4x4 and 8x8 matrices provide the same re- 
sults? This is due to the fact that the eigenvalues of both 
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FIG. 2. The direct counting results of ID for (a)R4, (b)R22 
and (c)R110 with different L. 
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FIG. 3. The In-degree Distributions for (a) R193, R147, 
R110 and (b) R137, R124 and R54 with L = 25. 
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FIG. 4. The In-degree Distributions for R110, R54, R73 and 
R62 with L = 25. 
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FIG. 5. The in-degree distributions for various rules with (a) 
v = 0.551 and (b) v — 0.6942; (c) The in-degree distributions 
for Class II, where v = 0.811 for R4, v = 0.3471 for R50 and 
v = 0.4650 for R218. All of results are ploted with L = 25. 



obtained in Amer Shreim et al.'s study [3[ is: 



matrices can be very much different and there is no obvi- 
ous reason for them to be the same. For all the cases we 
have studied the maximum in-degree can be exactly cal- 
culated by both 4x4 and 8x8 matrices and the results 
are identical. Unfortunately there exists no analytical 
proof. To illustrate this fact, the following example is 
the calculation of k m with the T matrices for Ris- Fig. 
[7] is the evolution rule of Ris, and the 4 x 4 T matrix 



/ 1 
11 
10 

Vo o i i 



(7) 



The hubstate of R18 is (00... 00), then Fig [7] shows 
that (001) and (100) are excluded from the pre-images 
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FIG. 6. The diagram for 8 x 8 T matrix. 
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FIG. 7. The diagram of R18. 



of (00 . . . 00), therefore one has Tig: 



/10000000X 
00000000 
1 
1 1 



Tie 



(8) 



00000000 
1 1 
1 

yooooooii/ 



The in-degree of hubstate is given by k(S) — Tr(T L ), 
with L = 10 one has: 
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It is clear that the traces of these two matrices are the 
same and resulting the same value for the maximum in- 
degree k m . 

We have also obtained the 8 x 8 T matrix for more 
than 22 rules and also rules related to them by symme- 
try( by exchanging o 1 ). They are: 1, 2, 4, 5, 6, 7, 
18, 19, 22, 23, 33, 36, 40, 54, 62, 73, 94, 110, 160, 218 
and 232. Except for R94, all these rules possess both T 
and T which results with the same k m . Further appli- 
cation of this augmentation procedure is illustrated by 



the example of the majority rule with more neighbor- 
ing cells included in the evolution. The evolution rule is 



R(4 



i s-s: 



s' +1 where 



q\ >0 
qj <0 

it 



= 



(11) 



and q\ — 53u=-2 S i+ M - The T matrix for this rule is a 

32 x 32 matrix and the maximum in-deg 1414 

for L = 15. From this example is is also apparent that, 

in general for the evolution rule involves 2n neighboring 

cells, the basic T matrix should be defined by 2 2n+1 x 
2 2n+i matrix . 

From the above discussion, one can see that with 8x8 
T matrix, more rules can be treated within this approach. 
As a consequence, such extension of T into larger matrix 
seems to provide a better description for network anal- 
ysis. Unfortunately such approach can still not be able 
to cover all elementary CA. There exists rules in which 
the configurations of hubstates vary according to L being 
even or odd, such as R50 and R77. On the other hand, 
there are also cases in which T varies with L indefinitely, 
one of the good examples is R43. We had analyzed R43 
with different L and the result is listed in Table HTT1 The 
number of hubstates will change with L, and the period 
of the attractors is large, therefore it is not appropriate 
to calculate k m by using T matrix. For these situations, 
there is no scale of k m with arbritrary L. 

For rules with the above irregularity special treatment 
is required for each particular case. In the following we 
will discuss R50 and R77 for illustration of their com- 
plication. The maximum in-degree of R50 and R77 can 
be properly described by 8 x 8 T matrices only when L 
is even, since the in-degree distributions will not depend 
with lattice length L. 

Fig. [8] shows the evolution according to R50, and Ta- 
ble IIVI shows the relation between maximum in-degree 
k m and lattice length L. The hubstate of R50 is cither 
(0101 . . . 0101) or (1010 . . . 1010) when L is even.The k m 

of R50 can be evaluated by k(SL=even) = l/2Tr[T 50 ], 
the factor 1 /2 is due to the fact that the 8 x 8 T matrix 
carries the information of all preimages of both hubstates, 
namely (0101 . . . 0101) and (1010 . . . 1010). For even a, 
all the preimages of L = a + 1 can be obtained by just 
appending "0" or "1" at the boundary to all the preim- 
ages of L = a. As a consequence, both even L = a and 
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FIG. 8. The diagram of R50. 



TABLE IV. Maximum in-degree corresponding to lattice 
length L in R50 by direct counting. 
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odd L = a + 1 have the same maximum in-degree k m . 
From these results one can see clearly that for even L k m 
follows scaling law with v ~ 0.345. But for odd L there 
is no scaling behavior for k m . The same conclusion also 
appears in R77 which is presentated in Table [V] 

Rule 77 is treated in a different way. Fig. [9] is the evo- 
lution according to R77. The network configurations are 
completely different for L being odd or even. For illustra- 
tion it is shown in Fig. [TU]the networks for L = 4 and 5. 
Hence one would expect k m varies with L. The hubstatcs 
of the R77 are (0101 . . . 0101) and (1010 . . . 1010) when 
L is even, but k m of odd L is quite varying, as given in 
Table El For R77, where L = 10 we observed that there 
exist only two hubstates where 22 hubstates are found 
for L = 11. Concentrating on one of L = 10 hubstates, it 
is clear that the hubstate of L = 11 is just adding 1 or 
to the end of L = 10's hubstate. Furthermore, there are 
totally 46 preimages for L = 10 and only 28 for L = 11. 
The 18 italic configurations with L = 10 do not associate 
with the preimages in L = 11. From Fig. [TT] it is known 
that only those preimages of L = 10 without (111) at the 
boundary for can associate to the preimages of L = 11. 

To evaluate the k m with odd L we follow the approach 
suggested in Q in obtaining arbitrary in-degree of R4, 
where apart from the T matrix itself an auxiliary matrix 
is constructed to resolve the difficulty. Therefore, by this 
approach we propose the following formula: 



k(S L = 



:>dd 



~L-1 

= Tr[T 77 



x T, 



multiplied] 



(12) 



This is due to the fact that the hubstates of odd L can be 
viewed as just adding "0" or "1" in the hubstates of even 
L. According to Fig. [TUl we could assign the T mu i tip i ied , 
which just corresponds to adding "0" in the hubstate 
obtained by referring to Fig. 1111 since the mapping at 
the boundary of preimages with L = 11 has only two 



State 








MJ3 




CM3 




cm 


Rule 77 


□ 


■ 


□ 


□ 


■ 


m 


□ 


■ 





FIG. 10. The networks of R77 with (a) L = 4 and (b) L = 5 
which are drawn by PAJEK. 



TABLE V. Maximum in-degree corresponding to lattice 
length L in R77 by direct counting. 



kjn 


even L 
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km 


odd L 
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46 


10 
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28 


11 


0.437032 
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13 


0.454376 
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15 


0.467415 
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16 
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17 


0.477279 
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Preimages of (10101010100) 
with L=ll 
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0101010101 
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FIG. 9. The diagram of R77. 



FIG. 11. The lists of pre-images of particular hubstates in 
R77 with different L. 
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TABLE VI. Summary for maximum in-degree calculation for 
R50 and R77. 



Rules 


Method 




R50 


k \Sl= even 


) = |Tr[Tg ] 






1 r~i-l, 




k(SL = odd) 


= |7>[T 60 ] 


R77 


k{SL = odd) 
k(S L = odd) 


= kTr[T L 77 ] 

r~L-l ~ 

— Tr [T77 x T mu iti p iied\ 



possibilities: 



^multiplied 

000 >= |000 > 

TmultipUed \ OH >= |H0 > 



(13a) 
(13b) 



Therefore we have T, 



multiplied 



as 



multiplied 



(14) 



/l \ 
00000000 
00000000 
1 
00000000 
00000000 
00000000 

\fl / 



Then we can get the exact maximum in-degree of hub- 
state with any lattice length L. This analytical results 
agree with direct counting for L < 25. A summary of the 
formulas for the maximum in-degree of R50 and R77 are 
given in Table IVI1 

In this section we have found more rules which could 
not be treated with 4 x 4 T matrix can now be dealt with 
by 8 x 8 transfer matrix. Consequently, the maximum in- 
degree and the ID can be evaluated efficiently. 



IV. CONCLUSION 

A detail discussion of the local characteristic of net- 
work generated by elementary CA is presented in this 
work. First we have located the possible source of er- 
rors in claiming the multiscaling of such network in 
This is due to a misused of derivative on large difference. 
Secondly the in-degree distributions of many rules were 
calculated by direct counting method and the multiscal- 
ing characteristic does not appear. From the studies up 
to L = 25, it suggests that the in-degree distribution of 
CA is L independent regardless of its classes. Further- 
more all the possible largest eigenvalues X m for 4 x 4 
transfer matrix T are obtained and analysis on the cor- 
relation between v — log 2 X m and classification of CA is 
discussed. We have found that there might be more rules 
sharing the same in-degree distribution but correspond- 
ing to different classification according to S. Wolfram. It 
has been found that there are good correlation for class 
IV which possess v = 0.4065, however, apart from that, 
the ID of other classes do not share any universal pattern 
and also the value of v can not be used as an indicator for 
classification. Therefore the in-degree distribution is not 
a proper characteristic to classification which was also 
pointed out in Q . In section IIII1 the transfer matrix was 
extended to 8 x 8 T matrix. This augmentation in some 
sense respects the evolution of CA rule and is also more 
natural to reflect the characters of hub-states' pre-images. 
As a result the 8x8 matrix scheme was applied to sev- 
eral rules which do not possess a 4 x 4 T matrix. The 
32 x 32 T matrix has also applied to CA beyond neighbor 
evolution. The L dependence of the hubstates of some 
particular rules such as R50 and R77 are discussed in 
detail. The calculation of their k m are also treated in 
this work. 
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